Quantum interferometry with and without an external phase reference 
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We discuss the role of an external phase reference in quantum interferometry. We point out 
inconsistencies in the literature with regard to the use of the quantum Fisher information (QFI) in 
phase estimation interferometric schemes. We discuss the interferometric schemes with and without 
an external phase reference and show a proper way to use QFI in both situations. 
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Laws of quantum mechanics impose fundamental 
bounds on measurement precisions of basic phys- 
ical quantities such as position, momentum, energy, 
time, phase etc. Theses bounds follow from the 
structure of the theory itself which contrasts the situ- 
ation encountered in classical physics where meas- 
urement uncertainties are due to factors which in 
principle may be eliminated by improving the qual- 
ity of measurement procedures. One of the most im- 
portant measurement techniques where such bounds 
have been analyzed is optical interferometry [1]. 

In a generic interferometric measurement using a 
Mach-Zehnder setup and classical light sources the 
precision of estimating the relative phase delay (p in- 
side the interferometer is bounded by the so called 
standard quantum limit (SQL) 5<p > where 
N is average number of photon-counts. At the fun- 
damental quantum level, the bound is a result of 
an independent probabilistic behavior of individual 
photon propagating through the interferometer. 

Breaching the SQL requires the use of special non- 
classical states of light were photons can no longer 
be regarded as independent. One of the first pro- 
posals in this direction was the idea to mix coher- 
ent light with the squeezed vacuum at the input 
beam splitter of the Mach-Zehnder interferometer 
[2]. Thanks to the reduced vacuum fluctuations in 
one of the quadratures of the squeezed state it is pos- 
sible to improve the precision beyond the SQL. This 
observations prompted the search for more funda- 
mental bounds on achievable precision, which would 
be obeyed by all quantum states [3]. 

In general, looking for the optimal phase estima- 
tion protocols is difficult since one needs to optimize 
over the input state \ipin) that is fed into the inter- 
ferometer, the measurement {II„} that is performed 
at the output and the estimator ip(n) — a function 
that assigns a phase value to a given measurement 
outcome. One of the popular ways to obtain useful 
bounds in quantum metrology, without the need for 




Figure 1. An interferometric scheme with coherent and 
squeezed vacuum states interfered at a beam-splitter, 
with arbitrary quantum measurement potentially in- 
volving an additional reference beam. In general, the 
QFI bounds on the phase estimation precision depend 
on the way the interferometer phase delay is modeled: 
(i) phase shift only in the upper arm, (ii) phase shift 
distributed symmetrically, (iii) phase shifts defined with 
respect to an additional reference beam. 



cumbersome optimization, is to use the concept of 
the quantum Fisher information (QFI) [4]. 

The purpose of this paper is to give a proper in- 
terpretation to the bounds obtained via the QFI and 
point out conflicting approaches where seemingly 
equivalent physical models lead to different quant- 
itative statements. We show that the source of the 
problem lies in the use of quantum states of light 
which are coherent superpositions of different total 
photon number terms without properly taking into 
account the role of an external phase reference beam. 

Let a, b be the anihilation operators of respect- 
ively upper and lower input modes of the inter- 
ferometer. For the purpose of this paper we con- 
sider the input state of the form |-0i n ) = \r, a), 
where \a) is the coherent state, a\a) = a\a), while 
\r) = exp[|r*a 2 — ir^) 2 ] |0) is the squeezed va- 
cuum state with squeezing parameter r (see Fig. 1). 
After it has evolved through the beam-splitter with 
power transmission r, and experienced the relat- 
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ive phase shift inside the interferometer U v , the 
state becomes \i/j v ) = U v B T \ipi D } where B T = 
exp[— i asm(s/r)(<vb + aw)], U v = exp[—i(paJa]. In 
a standard Mach-Zehnder setup one interferes the 
two modes on another balanced beam-splitter and 
detects number of photon clicks n a , n b in the two 
output modes. In an idealized setup with no losses, 
perfect interferometer and 100% detection efficiency 
this leads to a phase dependent probability distribu- 
tion of clicks: 



p(n a ,n b \ip) = \(n a ,n b \B 1/2 U v B T \'ipu 



(1) 



Instead of looking for the best possible estimator 
of the phase, which in general is a hard task, one 
can invoke the Cramer-Rao bound [5] which states 
that for k repetitions of an experiment and any loc- 
ally unbiased estimator ip(n a ,n b ) the uncertainty of 
estimation is bounded from below by: 



Sip > 



1 



i 



p(n a ,n b \ip) 



dp(n a ,n b \ip) 



dip 

(2) 

where F is the Fisher information. Moreover, the 
bound can be saturated in the limit k — > oo, by 
making use of the maximum likelihood estimator. 

A priori it is not obvious that this type of meas- 
urement is the optimal way to extract phase inform- 
ation from the state The quantum Cramer- 
Rao bound [4] provides an answer to this problem 
and states that whatever the measurement chosen 
the following bound on the estimation uncertainty 
holds: 



Sip > 



1 



F Q = 4(<vgv;>-i<v;ivv>i 2 ) (3) 



where \ip' v ) = , and Fq is called the quantum 
Fisher information. Fq depends neither on the meas- 
urement nor on the estimator and it is solely a func- 
tion of the probing state, which makes it an easy to 
calculate quantity. Moreover, one can always find a 
measurement (that may depend on the true value ip) 
for which F — Fq. In what follows, we will drop k 
for simplicity and use notation where Sip = 1/ \JFq. 

A lot of work in quantum enhanced interferometry 
has been based on utilizing the Fq [6-8] . In a typical 
approach one maximizes Fq over a class of input 
states satisfying some constraint (e.g. total energy) 
and in this way finds the input states optimal for 
quantum interferometry. 

Let us investigate the consequences of this ap- 
proach in our setup. The input state jV'in) = \ a , r ) 



has the mean number of photons equal to: n = 
\a\ 2 + sinh 2 r. Having fixed fi we look for optimal 
a, r and the transmission coefficient t that maxim- 
ize Fq . If one follows this procedure rigorously, then 
the solution depends strongly on the way the phase 
shift between the beams is modeled inside the inter- 
ferometer. As a simple illustration of this counter- 
intuitive behavior, the relative phase shift ip may be 
modeled in e.g. two ways depicted in Fig. 1 as (i) and 

exp[— iipa' a] 



(ii). These two cases correspond = t^m\—imn^r 



and U<£ % ' — exp[— i^a^a + i^b^b] respectively. When 
plugged into Eq. (3) they yield: 

F { q ] =4T 2 |a| 2 +2(l-T) 2 sinh 2 (2r)+S (4) 



F. 



(a) 



(l-2r) 



±sinh 2 (2r)l (5) 



where # = 4r(l — r) (|a| 2 e 2r + sinh 2 r) , and in or- 
der to simplify the formulas we have put the relative 
phase between the input beams to be 7r/2 (r — |r|, 
a = i\a\), which is the optimal choice for this and 
all the examples presented in this paper. The formu- 
las are clearly different, which becomes evident when 
we set t = 1/2, r = 0, in which case Fq* = 2|a| 2 , 
F^ = \a\\ 

To understand what lays behind this discrepancy 
consider an even more exotic case of r = 1, r = 0. 
The coherent state is simply transmitted to the up- 
per arm so there is no interferometer at all, yet 
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a\ . To give a meaning to 
these ,,unphysical" results, notice that QFI simply 
depends on the change of the probe state under the 
variation of the parameter ip. Even if we send a co- 
herent state | a) to the upper arm alone, then under 
the phase shift ip it evolves to |ae 4V ), which differs 
from | a) and in principle may provide us with useful 
information on the value of phase ip. The physical 
content that is missing in this reasoning is that the 
phase information is only available once we have ac- 
cess to an additional reference beam with respect to 
which the phase shift ip is defined. In other words 
there is no such thing as an absolute phase shift — a 
seemingly obvious fact which has nevertheless signi- 
ficant implications for the problem considered, and 
has been treated in contradicting ways in the liter- 
ature. 

The whole problem revolves around quantum 
states that are coherent superpositions of differ- 
ent photon number states such as e.g. coherent 
or squeezed states. Take a coherent state ||a|e l61 ). 
Since all measurements in quantum optics rely ul- 
timately on photon counts, no measurable con- 
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sequences of these coherences may be observed unless 
these states are interfered (as e.g. in a homodyne 
measurement) with a reference beam with respect 
to which the phase 9 is defined. Otherwise, one 
is entitled to phase average the state without any 
observable consequences, i.e. replace ||a|e* e ) with 
p = J ^\\a\e 10 } (\a\e te \ which is an incoherent mix- 
ture of photon number states with Poissonian stat- 
istics [9]. 

Going back to our quantum interferometric setup, 
if we indeed consider just the two modes of the inter- 
ferometer and do not allow any additional reference 
beam, then as an input we should rather consider a 
phase averaged state of the form: 

P(r ' a) = / ^ V o V oMM V o^ V o 1 (6) 

where Vg = cxp(— iOx^x). Notice, that squeezed and 
coherent states are averaged over a common phase 
9, which reflects the fact there is a physical mean- 
ing in the relative phases between them. Calculation 
of QFI for p(r, a) — F^ p) — are more involved since 
the state is mixed and instead of Eq. (3) one needs 
to employ a general formula involving the concept of 
the symmetric logarithmic derivative [4]. The res- 
ulting Fq" 1 is different both from Fq* and Fq 1 ^ and 
does not depend on the choice of the phase shift 
generator — be it U$ or u!p l K All that matters is 
the relative phase between the arms of the interfero- 
meter. Fq^ achieves maximum for r = 1/2 in which 
case it takes a simple form: 

maxfW = F^ r=1/2 = \a\ 2 e 2r + sinh 2 r. (7) 

There is a great deal of confusion in the literat- 
ure since formulas Fq\ Fq^ are often used instead 

of Fq ^ without discussing the need of an additional 
reference beam [8, 10, 11]. Despite this, one some- 
times arrives at the correct result, since e.g. for 
r = 1/2, F [ Q i] = F^> and that is why the results 
in [8] are indeed correct. However, had one used a 
phase shift generator (ii) instead of (i) one would ar- 
rive at a different solution. Similar objections can 
be raised in the context of [10] where Fq> is used 
and [11] where one defines standard quantum limit 
as Sip 2 = (2|a| 2 )~ 1 instead of (|a| 2 ) -1 which is again 
due to the use of F(p instead of F^ p) . Making use 

of Fq\ Fq 1 ^ without mentioning the need of a ref- 
erence beam is misleading since it is not clear what 
experimental setup these quantities really refer to. 



Let us now consider a situation in which we indeed 
have an access to an additional reference beam - 
represented by the state |/3) in Fig.l — and want to 
properly analyze the quantum interferometric setup. 
If the reference beam is strong, we can treat it as a 
phase reference for the other two modes. Therefore 
we introduce two phase shifts ipi, </? 2 as in (iii) in 
Fig. 1 which are defined with respect to the refer- 
ence beam. In a sense, we now face a two-parameter 
estimation problem. The proper way to proceed is 
to employ a two parameter Cramer- Rao bound [4]: 

S > T-\ Fij = AM^di^d^) - (c^XV^V)) , 

(8) 

where Xjj, i = 1, 2 is the covariance matrix for para- 
meters tpi, ip2, F is the quantum Fisher informa- 
tion matrix (QFIM), \i>) = V^V^ 2 B T \ip- m ) is the 
probe state after sensing the phase shifts <pi, ip 2 and 
\ditp) = fj^-- If one is now interested in the bound 
on the uncertainty of ipi the proper formula reads: 

6<Pi > ^(F- l )u- (9) 

Note that in general (F~ 1 )u ^ (Fa) -1 . 

In quantum interferometry we are interested in 
the phase shift difference between the interferometer 
arms i.e. <p_ = ipi — ip2, so it is more convenient 
to write QFIM in basis ip± = <p\ ± ip 2 . Calculating 
QFIM in ± basis yields: 

* - \ (l-2r)i5 (1-2t) 2 0+S ) (W) 

where © = |a| 2 +sinh 2 (2r)/2, Sj = sinh 2 (2r)/2-|a| 2 , 
and finally the bound on estimation precision of 
P- can obtained easily via Eq. (9). The minimal 
uncertainty is obtained for r = \ in which case 

5(p- > \J (|a| 2 e 2r + sinh 2 r) 

It is interesting to note that this is the same result 
as the one obtained for the phase averaged state us- 
ing Fq^ from Eq. (7). This observation proves that 
in the setup considered (with r = \) there is no ad- 
vantage in using the reference beam when estimating 
the phase difference between the two arms of the in- 
terferometer. More generally, it can be shown that 
this is a feature of all path-symmetric pure states, 
i.e the states that are symmetric with respect to an 
exchange of the arms of the interferometer [12]. It 
is also worth mentioning that the optimal measure- 
ment in our setup when r = 1/2, and more generally 
whenever we deal with a pure path-symmetric state 
in the interferometer, is a standard photon count 
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Figure 2. Bounds on the phase estimation precision cal- 
culated using the QFI, in case of an ideal (black) and 
lossy {n = 0.8, gray) interferometer. Different curves cor- 
respond to QFI calculated using different models: Fq* 

(dotted), -Fq"' (dashed), F Q p) (solid). In case of a lossy 
interferometer the additional reference beam may im- 
prove the precision: (J-^ 1 ) (gray, solid, bottom) < 

(-Fq^) (gray, solid, top) , while for the ideal interfer- 
ometer these quantities coincide. 

measurement after the two modes are interfered on 
a balanced beam splitter [7, 13]. 

As a summary of the discussion, in Fig. 2 we plot 
in black the bounds on Sip obtained using different 
QFIs. The bounds are plotted function of the 
total number of photons n used, and parameters (r, 
a, r) are chosen to maximize the respective QFI. One 
can easily notice that the uncertainties calculated us- 
ing Fq' and Fq 1 " 1 are overly optimistic. The reason 
behind this is an implicit assumption that e.g. in 
the case of Fq , the lower arm of the interferometer 
(where there is no phase shift element) is perfectly 
aligned with the reference beam. Such an assump- 
tion can hardly be justified in practice. 

Things become more complicated when one takes 
into account loss in the interferometer. Let rj be 
the power transmission coefficient in both arms of 
the interferometer. All results presented in the pa- 
per may be rederived in this setup although calcula- 
tions are more involved. Fig. 2 depicts in gray the 
resulting uncertainties for exemplary loss coefficient 
1 — r] = 0.2. Apart from a similar observation that 
Fq 1 and Fq 1 ^ yield overoptimistic results, we addi- 
tionally observe that (J 7-1 )-- < (F^ p) ) which is 

illustrated by a thick band and proves that having 
an additional reference beam helps in estimating the 



phase difference in a lossy interferometer. 

It is interesting to understand deeper what we 
really mean by strong reference beam. Clearly if 
|/3 1 is not strong enough we can hardly treat it as 
a phase reference. To solve this problem consider a 
phase averaged three mode state: 

-V e a V e b V£\r, a, 0)(r, a, (3\ V*V?V? . 

(11) 

Calculating the QFIM in this case can be done only 
numerically. Finally we can calculate the optimal es- 
timation strategy (optimal r, a, r) and the resulting 
bound on precision 6ip- > (J 7-1 )-- function 
of \(3\. With the increasing value of \(3\ we will ap- 
proach the regime discussed before where we treated 
the reference beam as strong enough so it can serve 
as a perfect phase reference. In the case of an ex- 
ample depicted in gray in Fig. 2 this corresponds to 
improving the estimation precision by going form the 
upper to the lower boundary of the gray band with 
increasing |/3|. 

A deeper analysis [12] shows that a sufficient con- 
dition for treating the reference beam as a perfect 
phase reference is |/3| 2 » n 2 . The reference beam 
needs to have much more than the square of the num- 
ber of photons traveling in the proper modes of the 
interferometer, a fact observed also in [14]. 

In summary we have pointed out some possible 
flaws in the interpretations of the results obtained 
using the QFI for states which are superpositions 
of different total photon number terms and showed 
that the full understanding of the problem is only 
possible if the role of an additional reference beam 
is properly taken into account. 
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